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Coastal  Engineering 
Technical  Note 


IRREGULAR  WAVE  ORBITAL  VELOCITIES  AT  COASTAL  STRUCTURES 


Purpose:  The  purpose  of  this  note  is  to  describe  an  available  PC  program  for 
estimating  near-bottom  water  particle  velocities  which  contribute  to  the 
problem  of  scour  at  and  around  various  coastal  structures  (including  vertical 
walls,  sloping  revetments,  and  rubble-mound  structures). 


Background :  Undermining  of  the  toes  of  coastal  structures  which  may  occur  as 
a  result  of  scour  hole  development  is  a  serious  problem.  The  undermining 
causes  damage  to  structures  and  often  to  the  buildings/beaches/roads  they  were 
designed  to  protect.  Coastal  engineers  have  long  recognized  the  consequences 
of  scour  at  the  base  of  these  structures,  and  expensive  toe  protection  schemes 
often  have  been  implemented.  In  instances  where  appreciable  scour  already  has 
occurred,  a  common  solution  has  been  to  fill  the  scour  hole  with  stone  or 
other  suitable  material.  Scour  potential  at  the  toe  of  coastal  structures  can 
be  assessed  if  the  magnitude  of  the  wave  orbital  velocities  can  be  estimated. 
An  estimate  of  irregular  wave  near  bottom  velocities  is  possible  using  a~new 
technique  developed  at  CERC.  Empirical  relationships  for  the  reflection  - 
coefficient  and  the  reflection  phase  that  result  from  irregular  waves  being 
partially  reflected  by  coastal  structures  can  be  used  In  conjunction  with  a 
theoretical  expression  (Hughes  1992)  to  estimate  the  horizontal  and  vertical 
root -mean- squared,  rms,  water  velocities  anywhere  in  the  water  column  seaward 
of  the  structure  (particularly  near  the  bottom  in  the  vicinity  of  the  toe) . 

The  theoretical  formulations  are  given  below: 
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The  above  equations  represent  the  approximation  of  the  water  particle 
velocities  under  random  waves  by  the  linear  superposition  (summation)  of  many 
linear  components,  where  the  subscript  n  refers  to  the  nth  individual 
component .  Also , 
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root-mean-squared  horizontal  water  velocity 

root -mean- squared  vertical  water  velocity 

vertical  coordinate  with  z- 0  at  the  SWL  and  z  -  -h  at  the 

bottom 

horizontal  coordinate  with  x-0  at  structure  toe 
water  depth 

gravitational  acceleration 

wave  number  of  the  nth  linear  component 

angular  wave  frequency  of  the  nth  linear  component 

incident  amplitude  of  the  nth  linear  component 

reflection  coefficient  of  the  nth  linear  component 

phase  shift  of  the  nth  linear  component  due  to 

reflection 


These  equations  were  developed  assuming  nonbreaking  irregular  waves 
propagating  over  a  flat  bottom  normally  incident  to  the  reflective  structure. 
The  rms -velocities  can  be  estimated  at  any  position  (x,z)  given  the  incident 
wave  spectrum,  water  depth,  and  the  reflection  and  phase  functions. 


Reflection  Phase  Relationship 


Analysis  of  flume  test  results  using  monochromatic  waves  imply  a 
relationship  for  phase  in  terms  of  an  Irr£barren-type  parameter,  namely 


^  ■  H 
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where  T„  is  the  wave  period  of  the  nth  component  and  tan  fl  is  the  slope  of  the 
seaward  face  of  the  structure  ( p  is  the  angle  the  slope  makes  relative  to 
horizontal) .  The  empirical  (best- fit)  equation  derived  from  this  analysis  is 
given  for  both  smooth  and  rubble-mound  structures  by 

0om 4*  (l  -  2.20^15)  (4) 

where  has  units  of  radians,  as  shown  in  Figure  1.  This  particular  form  of 
the  phase  equation  was  chosen  so  that  6n  *♦  4*  as  £h  -♦  0,  one  of  the  limiting 
cases  for  complete  reflection.  (Note:  Irregular  waves  cannot  be  used  to 
determine  the  reflection  phase  relationship  due  to  the  2k  cyclic  nature  of  the 
phases . ) 

REFLECTION  PHASE 
Sloping  ttuoteM 


Figure  1.  Reflection  phase  relationship  for  smooth  and  rubble  mound 
structures  for  non-breaking  waves. 
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Reflection  Coefficient  Relationship 

Estimates  of  rms -velocity  are  for  irregular  waves,  so  data  from 
irregular  wave  tests  were  used  to  establish  the  reflection  coefficient 
relationship.  Analysis  of  data  obtained  from  irregular  wave  conditions 
covering  3  water  depths,  3  peak  wave  periods,  smooth  plywood  slopes  of  1:2, 
1:3,  and  1:4,  and  a  rubble -mound  slope  of  1:2,  produced  the  relationships 
shown  in  Figures  2  and  3.  Once  again,  the  Irribarren-type  parameter,  , 
proved  to  be  a  good  nondimens ional  parameter.  The  data  shown  on  Figure  2 
yielded  a  best- fit  equation  given  by 


K  *  -i— P'_ —  ( Smooth  Structures) 
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which  is  shown  as  the  solid  line  on  Figure  2,  whereas  the  data  in  Figure  3 
produced  the  following  relationship: 
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Figure  2.  Reflection  coefficient  relationship  for  smooth  plane  sloping 
structures . 


REFLECTION  COBTOENT 


Figure  3.  Reflection  coefficient  relationship  for  rubble-mound 
structures . 
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The  above  empirical  expressions  for  reflection  coefficient  (Eqns.  5  and  6)  and 
reflection  phase  (Eqn.  4)  have  been  programmed  into  a  PC  SHAREWARE  program 
that  solves  Eqns.  1  and  2  for  given  input  incident  wave  spectra. 

Using  this  SHAREWARE  program,  estimates  were  made  of  the  horizontal  and 
vertical  rms -velocities  at  mid -depth  using  the  resolved  incident  spectra  from 
laboratory  tests  mentioned  above.  These  estimates  were  compared  to  actual 
measurements  of  irregular  wave  rms -velocities  as  shown  for  two  cases  in  Figure 
4.  The  structure  toe  is  situated  at  position  x-0,  and  this  represents  the 
closest  distance  to  the  structure  for  which  estimates  can  be  made.  Note  that 
the  location  of  the  maximum  u_ >  velocity  occurs  at  different  positions 
relative  to  the  structure  toe  tor  the  two  different  wave  conditions.  This 
illustrates  the  importance  of  including  the  phase  relationship  in  estimates  of 
potential  scour  at  the  toes  of  structures. 

The  following  sections  provide  details  about  the  PC -computer  SHAREWARE 
program  that  allows  the  user  to  estimate  the  rms  horizontal  and  vertical 
velocities  using  the  above  described  methodology.  Included  is  the  procedure 
for  using  the  SHAREWARE  and  generating  program  output. 

PC  Requirements 

The  present  version  (Ver  2.00)  of  the  code  has  been  compiled  to  execute 
on  a  PC  with  a  math  coprocessor.  It  can  be  used  immediately  on  systems  which 
have  a  math  coprocessor,  VGA  graphics  card,  and  a  HP  LaserJet  III  installed  on 
the  LPT1  port  (only  necessary  to  get  plot  hardcopies) .  Different  PC  setups 
will  need  to  install  the  program.  Instructions  for  installing  the  program  are 
very  straightforward  and  are  included  with  the  software. 

Executing  the  Program 

The  program  is  started  by  typing  RXSVEL  <return>  .  The  program  first 
examines  the  current  directory  for  a  default  file  named  RMSVEL.CFG  which 
contains  the  parameter  defaults  from  the  last  time  the  program  was  run.  If 
this  file  is  found,  the  default  parameters  are  read  into  the  program.  If  the 
file  does  not  exist,  hard-wired  defaults  are  used.  At  the  end  of  the  program, 
the  file  RMSVEL.CFG  is  created  and  user-modified  defaults  are  written  into  the 
file  for  use  the  next  time.  This  feature  saves  time  during  repeated  applica¬ 
tion  of  the  program  where  only  one  or  two  of  the  parameters  are  altered. 

User  Input  Section 


The  first  screen  (Figure  5)  is 
typical  of  the  interactive  screens 
which  prompt  the  user  for  input  or 
decisions  required  to  run  the 
program1.  This  screen  lists  the 
assumptions  inherent  in  the  analysis 
performed  by  the  program.  The  next 
screen  presents  the  first  option  for 
the  user  -  which  type  of  irregular 
wave  condition  to  use  as  input. 

There  are  two  basic  options  for  the 
input  incident  wave  spectrum: 


This;  program  calculates  the  uaue-irfueed  roo'-npnn-sqjared 
I  In id  velocities  seauard  cr  thp  toe  cT  a  coastal  structure 
urder  the  following  assumption:;: 


<1)  Hopnc.l  untie  incident  ( tuo-f!  irensi:>n j  1 ) 


12)  Honhrpakirn  irregular  nave? 

(3)  Horizontal  sea  fleer  ceauard  of  the  structure 


tob  choic: 
ENTER  SELECT 


Figure  5.  Initial  interactive  screen 
displayed  in  Shareware  program. 


1  Note:  Help  screens  are  normally  available  by  pressing  the  FI  key. 

4 


RMS  Velocities  -  irregular  Wave  Condition  1 
Smooth  i-on-2  Slope  (Dio) 


i  • 

•2  7 


-40  -20 


— UnM  Msaxur«d 

-  Unm  Calcutatwd 

—  wrnt  Measured 
-  Wnw  Calculated 


20  40  60  80  100  120  140  160 

Distance  From  Slop*  To*  •  (om) 


RMS  Velocities  -  Irregular  Wave  Condition  2 
Smooth  1-on-2  Slope  p14) 


-40  -20 


»  Unm  Measured 

-  unm  Caleulatad 

— wim  Maaaurad 
- Wim  Caleulatad 


20  40  60  80  1  00  120  140  1  60 

Distance  From  Slope  Toe  -  (cm) 


Figure  4.  Comparison  of  predicted  values  to  measured  values  of 
v,,,,  for  two  different  conditions  measured  in  the  laboratory. 


The  TMA  option  will  probably  see  the  most  use  by  the  field  offices.  This 
option  generates  an  incident  TMA.  spectrum,  which  is  a  spectral  representation 
of  waves  which  field  experience  has  shown  to  exist  with  specified  conditions 
of  depth,  peak  frequency,  and  windspeed.  The  TMA  spectrum  is  generated 
internally  by  the  program,  using  the  user -specified  parameters  H«,,  Tp,  and  7. 
Hno  is  known  as  the  zero  moment  wave  height,  and  is  the  energy-based 
definition  of  significant  wave  height,  H,.  In  deep  water,  1L,  is 
approximately  equivalent  to  H.,  but  the  two  can  be  quite  different  as  waves 
shoal  (Hughes  and  Fowler,  1990).  T_  is  the  peak  spectral  period  and  7  is  a 


spectral  peakedness  parameter  typically  ranging  between  1  and  7.  The  default 
7  is  3.3,  which  will  generate  a  spectrum  typical  of  that  associated  with  the 
widely  known  Joint  North  Sea  Wave  Program  (JONSWAP)  data.  The  user  also 
supplies  a  name  for  the  output  file. 

2)  User- Supplied  Spectrum  Option. 

The  user -specified  spectrum  option  allows  input  of  a  measured  (or 
otherwise  specified)  incident  spectrum  for  use  in  the  calculation.  The 
following  options  are  presently  supported: 

Generic:  Input  file  contains  values  of  frequency  (Hz)  and  continuous 
energy  density  (a  -s). 

RAF:  These  are  special  files  created  using  CERC's  laboratory  reflection 
analysis  procedure  on  Laser  Doppler  Velocimeter  (LDV)  data.  They  are 
obtained  as  output  from  the  program  REFLECT . 

GAF:  These  are  also  special  files  created  using  CERC's  PC-based  GODA 
analysis  program  PCGODA. 

For  each  of  these  cases,  reading  of  the  input  files  continues  until  an  End- of - 
File,  <EOF>,  mark  is  encountered. 

All  of  the  above  options  create  an  output  file  having  the  same  name  as 
the  incident  spectrum  input  file,  but  having  an  ".RMS"  extension.  For  program 
evaluation  purposes,  an  input  file  named  DIO. RAF  is  provided  with  the 
shareware  program  for  use  with  the  second  option. 

Existing  Output  Filenames 

If  the  program  detects  that  the  selected  (TMA  option)  or  generated 
(User- selected  option)  output  filename  already  exists  in  the  current 
directory,  a  screen  will  appear  which  offers  following  choices: 

Overwrite  the  existing  file 
Rename  the  existing  file 
Enter  a  new  output  filename 
Exit  the  program 

Input  Parameter  Specification 

The  final  set  of  options  are  the  same  irrespective  of  what  type  of 
incident  spectrum  is  chosen  to  produce  the  velocities.  A  menu  is  offered 
giving  the  user  options  to  select: 

Water  Depth:  Depth  at  the  toe  of  the  structure  (in  meters). 

Solution  Depth:  Depth  at  which  velocities  are  calculated  (in  meters).  This 
entry  should  be  the  distance  down  from  the  SWL  in  meters.  For  example,  the 
bottom  would  be  the  same  as  the  depth.  It  is  necessary  to  include  a  negative. 
Also,  attempts  to  enter  a  value  greater  than  the  depth  will  not  be  accepted. 

Horizontal  Increment:  The  program  first  calculates  the  wavelength  associated 
with  the  peak  period,  then  divides  the  distance  of  three  of  these  wavelengths 
into  equal  length  increments  at  which  solutions  will  be  computed.  The  default 
for  the  number  of  horizontal  increments  is  200,  but  the  user  may  input  any 
value  within  the  range  of  1  -  200.  Units  for  horizontal  distance  are  meters. 

Structure  Slone:  Structure  slopes  are  entered  as  frsctions  representing 
rise/run.  Because  the  mildest  slope  represented  in  the  experimental  data  was 
1/4  (slope  -  0.25),  acceptable  slope  values  range  between  0.25  <  tan  <  •  . 
The  <o  condition  represents  case  of  a  vertical  wall. 
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Structure  Type:  Structure  type  can  be  either  smooth  slope  or  rubble-mound 

slope.  The  empirical  relationship  for  rubble-mound  slopes  is  based  on  a  slope 
of  1:2,  but  is  thought  that  reasonable  estimates  will  be  produced  for  other 
typical  rubble -mound  slopes. 

NOTE:  To  input  the  incident  spectrum  contained  in  file  DIO. RAF,  the  following 
parameters  are  associated  with  this  wave  condition: 

Depth:  h  —  0.305  m 
Slope:  tan  -  0.5 

Calculation  Section 

After  receiving  correct  user  input,  the  program  first  calculates  the 
continuous  incident  wave  spectrum,  then  offers  the  user  the  option  of  plotting 
resul ts  on  the  screen.  The  user  can  zoom  the  plot,  or  print  the  plot  if  the 
program  is  correctly  installed  for  the  printer.  An  example  of  a  plot 
calculated  for  a  TMA  spectrum  generated  using  H^,  -  1.20  m,  Tp  -  10.0  sec, 
and  7  »  3.3,  is  shown  below  in  Figure  6. 


INCIDENT  UAUE  CONTINUOUS  TMA  SPECTRUM 


Figure  6.  Example  of  TMA  Spectrum. 

The  next  step  is  the  calculation  of  the  rms  -velocities  at  each  of  the 
200  horizontal  locations.  A  counter  on  the  screen  lets  you  know  how  the 
computation  is  progressing.  Once  the  computing  is  completed,  the  spatial 
variation  of  rms -velocities  as  a  function  of  distance  from  the  toe  of  the 
structure  can  be  viewed  on  the  screen  and/or  printed.  Examples  of 
rms -velocity  plots  associated  with  the  inputted  TMA  spectrum  are  shown  in 
Figures  7  and  8.  Figure  6  presents  results  for  a  mid-depth  location  and 
Figure  8  presents  calculated  velocities  at  the  bottom.  Note  that  vertical 
velocity  is  zero  near  the  bottom. 
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ROOT-MEAN  SQUARED  IRREGULAR  HAVE  VELOCITIES 

Output  File:  TNA6FZ.RHS 

Swooth  Structure  with  Slope  *  6.500 

RTISUEL  Progrtn  -  Version  2.00 


Figure  7.  Plot  showing  calculated  velocities  for  mid-depth  condition. 


ROOT-MEAN  SQUARED  IRREGULAR  UAUE  VELOCITIES 

Output  File:  HM0FZ.RNS 

Swooth  Structure  with  Slope  *  6.566 

RNSUEL  Progrew  -  Version  Z.66 


Figure  8.  Plot  showing  calculated  velocities  for  near  bottom  location. 
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Output  Section 

After  leaving  the  plot  option,  a  final  screen  is  displayed  showing  the 
parameters  for  this  calculation.  This  same  information  is  written  to  the 
output  file,  along  with  the  200  values  of  horizontal  and  vertical 
rms -velocity .  A  truncated  sample  of  the  ASCII  output  file  is  shown  below. 


ESTIMATION  OF  ROOT-MEAN-SQUARED  IRREGULAR  WAVE  VELOCITIES 
■  ■  tit >  Program-Generated  TMA  Spectrum  with  GAMMA  *  3.30 

Version  1.00  Analysis  Performed  on:  07/07/93  at  10:26:14 


>  USER-SPECIFIED  PARAMETERS: 

Incident  Significant  Wave  Height  (Hmo) . . . . [  1.2000]  m 

Incident  Peak  Special  Wave  Period  (Tp) . . .  [  10.0000]  s 

Water  Depth  (h) . . . - . [  5.0000]  m 

Velocity  Solution  Depth  (z) . [  -4.9000]  m 

Structure  Slope  (rise/run) . [  0.5000] 

Structure  Type. Smooth  ] 

>  MAXIMUM  RMS-VELOCITIES 

Maximum  Urms  «  0.6133  m/s  at  Location,  X  «  0.0000 
Maximum  Wrms  *=  0.0068  m/s  at  Location,  X  *  16.2398 


X  Urms  Wrms 

(m)  (m/s)  (m/s) 


0.000 

1.015 

2.030 

3.045 

4.060 

5.075 

6.090 

7.105 

8.120 

9.135 

10.150 

11.165 

12.180 

13.195 

14.210 

15.225 

16.240 

17.255 

18.270 

19.285 

20.300 

21.315 

22.330 

23.345 

24.360 

25.375 


0.6133 

0.6032 

0.5903 

0.5749 

0.5571 

0.5375 

0.5165 

0.4950 

0.4736 

0.4532 

0.4344 

0.4177 

0.4035 

0.3922 

0.3842 

0.3799 

0.3796 

0.3832 

0.3907 

0.4016 

0.4150 

0.4302 

0.4464 

0.4628 

0.4788 

0.4941 


0.0055 

0.0056 

0.0057 

0.0058 

0.0060 

0.0061 

0.0063 

0.0065 

0.0066 

0.0067 

0.0067 

0.0067 

0.0067 

0.0068 

0.0068 

0.0068 

0.0068 

0.0068 

0.0068 

0.0067 

0.0066 

0.0065 

0.0064 

0.0063 

0.0062 

0.0061 


Near-Bottom 
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SUMMARY:  Engineering  approximations  for  the  mis  horizontal  and  vertical  velocities  in  front  of  coastal  structures 
can  be  estimated  in  terms  of  structure  slope,  water  depth,  and  incident  wave  spectrum.  The  formulations  in  this 
technical  note  are  based  on  linear  wave  theory  and  should  provide  the  practicing  coastal  engineer  with  useful  tools 
for  calculating  both  vertical  and  horizontal  velocities  under  irregular  waves  at  a  coastal  structure.  These  formulations 
have  been  verified  in  the  laboratory  but  have  not  yet  been  verified  with  field  data.  However,  the  estimates  they 
provide  should  be  sufficiently  accurate  so  long  as  wave  nonlinearities  are  minor.  For  nonlinear  wave  conditions 
(breaking  waves),  these  estimates  must  be  considered  preliminary. 

ADDITIONAL  INFORMATION:  For  additional  information  contact  Dr.  Jimmy  E.  Fowler  of  the  Coastal  and 

Hydraulics  Laboratory. 
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